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Abstract 

The tangential layers are characterized by a bulk plasma velocity 
and a magnetic field that are perpendicular to the gradient direction. 
They have been extensively described in the frame of the Magneto- 
Hydro-Dynamic (MHD) theory. But the MHD theory does not look 
inside the transition region if the transition has a size of a few ion gy- 
roradii. A series of kinetic tangential equilibria, valid for a collisionless 
plasma is presented. These equilibria are exact analytical solutions of 
the Maxwell- Vlasov equations. The particle distribution functions are 
sums of an infinite number of elementary functions parametrized by a 
vector potential. Examples of equilibria relevant to space plasmas are 
shown. A model for the deep and sharp density depletions observed 
in the auroral zone of the Earth is proposed. Tangential equilibria are 
also relevant for the study of planetary environments and of remote 
astrophysical plasmas. 
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I. INTRODUCTION 



The Swedish Viking spacecraft orbiting in the auroral region of the magne- 
tosphere discovered deep plasma density depletions [Tj. Their size along the 
magnetic field direction reach thousands of kilometers. They can be con- 
sidered as channels parallel to the magnetic field. These observation were 
confirmed by measurements onboard the Fast Auroral SnapshoT explorer 
(FAST) spacecraft [2]. The data showed a lot of plasma activity in and 
around these regions, particularly plasma acceleration (at the origin of the 
polar auroras) j3j, electromagnetic radiation [3] and a high level of electro- 
static turbulence. Most of the studies focused on the acceleration and 
the nature of the turbulence [7j, |H] 0- But surprisingly, little work has been 
devoted to the structure of these plasma cavities. Some theoretical works 
have been devoted to their origin, but none has tackled the question : are 
auroral plasma cavities in equilibrium, or do they vanish a soon as the cause 
of their formation disapears ? 

The aim of the work presented in this paper is to construct a kinetic 
model of the auroral plasma cavities. But as a very general solution to this 
problem is found, the field of applications appears much wider. I present 
here a family of equilibrium solutions that can be used in many domains of 
the physics of plasmas. They will be named as tangential equilibria because, 
as in the MHD tangential discontinuities, the magnetic field and the bulk 
velocities are tangential to the plane (y, z) of invariance. 

Previous works on tangential equilibria JU] , jH| have shown the existence 
of simple isothermal solutions that can be computed analytically. Other 
works, mainly focused on the study of the Earth magnetopause [12], have 
shown solutions that satisfy a larger class of constraints, but they cannot be 
constructed analytically. 

Most of these equilibria a based, for each species, on distribution functions 
that are products of an exponential function of the total energy and of p 2 z , 
by an arbitrary g function of p y , 



where a z = m/2T z = l/vf z , a± = m/2T± = l/v^ ± , are the reciprocals of the 
squared thermal velocities ; v z = p z /m and v\ = v 1 + Vy. If g(oo) ~ 1, then 
/ is a bi-maxellian distribution for x ~ oo. 




(1) 
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Krall and Rosenbluth, have derived a solution of the linearized ki- 
netic equations, by setting g(p y ) as a linear function. They used this simple 
solution for the study of gradient instabilities. 

Many authors (reviewed by Roth et al. [121), w h° studied tangential 
discontinuities, considered solutions with g functions that involve error func- 
tions. Most of these studies where developed for the study of the equilibrium 
of the Earth magnetopause. They include an electric field E x and no charge 
neutrality. There is no analytical solution. The differential equations on E x 
and A are solved through a numerical integration. 

Channell ^T] considered solutions with no electric field. He showed that 

A more general formula will be shown in Sec. II. 

In the same paper, Channell studied a few examples of solutions with 
simple but nonlinear g functions. In particular, he showed that there exists 
a solution with g(p y ) = exp (—rjpy), where 77 is a constant number. 

Harris jTU] built an equilibrium model of a current sheet. His model was 
widely used later for the study of the tearing mode instability. The Harris 
current sheet equilibrium corresponds to the case g{p y ) = exp iyp y ). 

All the solutions mentionned above are based on elementary g functions. 
They do not let the possibility of setting a priori a density profile, as in the 
more simple MHD and bi-fluid models. We can build a distribution function 
that is the sum of such elementary distribution functions and assert that 
each elementary part of the distribution corresponds to a family of trapped 
particles. 

Because of the nonlinearity of the equations, the density and temperature 
profiles associated to a sum of families of trapped particles are not simple 
combinations of the profiles associated to single elementary particle distri- 
butions. Nethertheless, using a numerical equation solver, Roth et al. [12 
added up to four elementary distributions functions (all of them were built 
with error functions). Differential equations were solved numerically. 

I show in this paper how to combine an infinite number of elementary 
solutions and build a plasma that is made of a continuum of famillies of 
trapped particles. The plasma equilibria presented in this paper are com- 
puted analytically. They do not have the high degree of generality of those 
shown in Roth et al. [12], (because I impose a null electric field), but they 
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are more general than those of Harris ^U] > and Channell ^T] . In particular, 
the plasma does not have to be isothermal. 

The equations of the kinetic model are developed and solved in Sec. II. 
The properties of the solutions and some particular cases are briefly analysed 
in Sec. III. An equilibrium model for the deep non isothermal density deple- 
tions of the Earth auroral zone is presented. The last section is a discussion 
about other applications and further developments. 

II. FORMULATION OF THE MAXWELL- 
VLASOV MODEL 

We consider a monodimensional equilibrium : dt — and V = (d x ,0, 0). 
We choose the z direction along the constant direction of the magnetic field 
B = (0, 0, B z (x)), (V • B — is therefore ensured). The magnetic field 
derives from a vector potential A(x) = (0, A y (x),0) such that B z = d x A y . 
This vector potential satisfies the Lorentz gauge (that is V-A = in the 
stationary case). The x direction is called the normal direction, the y, z 
plane is the tangential plane. 

In a monodimensional equilibrium, the particles are confined between two 
points where v x — 0. In the case of a null electric field, the size Ax of this 
area is 

Ax = 2v±f <oo c > (3) 

where < lu c >= (q/mAx) J5 2 B z (x). Absolutely no restriction on the varia- 
tions of B x is made. 

As we look for a monodimensional equilibrium, the total energy E = 

2 

+ q& of any particle and the generalized momentum p y = mv y + qA y and 
p z = mv z are the invariants of the motion. Any distribution function of the 
form / = f(E,p y ,p z ) is a solution of the Valsov equation. All the kinetic 
solutions studied in this paper will be of that form. The other character 
imposed on the solutions is a null electric field, and therefore no charge 
separation. 

The basic idea is to decompose linearly the g function (already men- 
tionned in Eq. ^) over a set of elementary distribution functions that corre- 
spond to an analytical equilibrium solution. As will be shown in this section, 
if the sum is parametrized by a parameter a that correspond to shifts of the 
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elementary distribution functions in the space of the potential vectors (and 
not in the space of the configurations), there is an analytical solution. 

For instance, we can show that a continuous linear combination of Chan- 
nell and Harris-like g distribution functions given by 

g(py) = n + / da n„(a) exp [—ri(a)(— a) 2 + v(a)(— a)], (4) 

J ai mm mm 

leads to an equilibrium that can be computed analytically. 

The functions n g (a), rj(a) > 0, and u(a) are defined almost arbitrarily 
and n is a constant scalar. Each of these three functions provide a degree 
of freedom and it is possible, as in MHD and bi-fluid models, to control the 
plasma density profile. But the distribution of the particles depends on the 
energy through an exponential factor that does not give any freedom to set 
the temperature profile. Therefore, instead of using [T] and EJ let us set : 

/ = P da ( Q « ^ ) l/2 e (-a, («)(„, -u z )*- a± (a K ) (wp + (5) 

and 

g a (p y ) = n g (a)e~^^-^ 2+ ^^-^. (6) 

With these distribution functions, a control of the plasma density profile 
(through n g (a)) and of the temperature profile of each species (through a z (a) 
and a±(a)) is possible. The parameter a is homogeneous to a vector potential. 

In order not to overload the equations, we do not express systematically 
the dependance of the parameters on a and x hereafter. 

To solve the charge neutrality equation (rij = n e ) and the Ampere equa- 
tion 

d 2 

~ j-t A v = VoJy, (7) 
we need to compute the contribution of each species s to the charge density, 

n s = da ( — ) ' 2 / dv y exp (— a±v 2 )g(mv y + qA y (x)). (8) 

Jai 7T J 

and to the current density, 

J ys = J da q(-^-) 1/2 J dv y v y exp (-a_LV 2 )g(mv y + qA y (x)) (9) 
daq(- ) 1 ^ 2 [ dv y exp (— a±v 2 )g'(mv y + qA y (x)) (10) 



"Ana 
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where g' is the derivative of g. These charge and current densities can be 
computed explicitely. We reorder Eq. (jHJ and (JBJ to separate the terms that 
depend on v y and v y , 



f = C da J^e-^*-^-^ 2 x (11) 

{n e~ a±v y + n (a)e~^ a ^™( A y&)- a )?+ v {<^ 

with 

1 Q 

P(a) = a ± + r]{a) et Q(a) = -u(a) - rj(a) — (A y (x) - a) (12) 

The contribution of each species s to the charge and to the current densities 
are 

n s = n 



+ ["'dan J^ e -vWrn)>At SopQ 
Jai V 7T 

where 5o,p,q is defined by 

So,p, Q = [ + °°dv y e^ +2 ^ = 

J — OO V 1 



and 



where 



Hence, 



J ys = r da qn g ^e-^ 2A yS 1 , P , Q , 

Jai V 7T 



Si, 



,P,Q 



^d W <-**+***> = |^ exp(^; 



n + da n„ 



J, 



a2 . n n 



_g 4(q:j_+?7) g L m 



-2-(A,-a)][r,£(^-a)-i,] 



(13) 

(14) 
(15) 
(16) 

(17) 



-g / da — 

Jen 2 y aj_ + 77 

q 



^ 1 

g4(a ± +rj) y 



[-277— (A, - a) + u]e -x+" 1 ™ 
m y 



a± + r] 

(A y -a)][ V %(A y -a)-v) 



Let 



e(«) 

5(a) 



a±{a)r}(a) q_ 2 



a±(a) + 77(a) m 
a_i_(a)u(a) . q 



a±(a) + 77(a) m' 
iVo(a) = n g (a), 



a±(a) 



\ a j. (a) + r/(a) 



g 4(a ± +?7) ^ 



(18) 

(19) 
(20) 

(21) 
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The contribution of each particle species to the particle density is: 



n s (x) = n + C da N e-^- a)2+S{Ay ~ a) (22) 



a i 



Let us set 
then 



na { X ) = No ( a ) e -tta)(Ay(x)-af+5(a)(A y (x)-a)^ ^ 



ra.2 

n s (x) = n + / dan a (x) (24) 

J a\ 

If the following equations are satisfied for each value of a, 

^('3')jons = £(^0 electrons (25) 
$(a)ions = 5(a) 

electrons 

(26) 

^o( a )ions = N (a) electrons, (27) 

then, the charge neutrality is satisfied. It is easy to show that one can freely 
choose the functions a±, r], v, n g that correspond to a particle species, then 
compute the £, 5 and iVo functions and find the 7], u, n g functions of the 
other species that are associated with the same £, 5 and iVo- There is only 
one restriction: rj must be positive. If we choose electron parameters, the 
condition for a positive ion rji parameter is, for each value of a 

'"' >2[(^)T ±t (a)+T ±e (a)} (28) 



r) e (a) m e 

which requires small enough values of i] e (a). Noticing that r\ homogeneous 
to the reciprocal of a squared velocity rj = v%, the above relation writes 

< > C£t + ^te (29) 

IIL e 1 e 

Let us call rj c (a) = v c (a)~ 2 the value that corresponds to an equality in Eq. 
()28|) and (|2^j). The case of r\ e = r\ c corresponds to an infinite value of rji. 
We shall see in the applications that because of the condition r\ < r) c , v v can 
reach very high values, and it is necessary to remember that v v is not an 
actual velocity, it does not correspond to a propagation phenomenon. 
Let us define 

k{a) = -fM)N (a)( e -- + -M = -2u iV (a)(T ±e (a) + T ±i (a)). (30) 

a_Le(a) a ±i {a) 
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Negative values of k(a) correspond to a positive value of n(a), and therefore 
imply an increase of the plasma density. On the contrary, positive values of 
k(a) induce a reduction of the density. The total current density is the sum 
of the ion and the electron currents, 

J y {x) = -— f 2 da ~[-2Z(Ay - a) + 5] e -£(A,(s)-«) 2 +W*)-«). (31) 

fJ,Q J ax 2 

Let us notice that Eq. If22j) . (j^ . and (jSH) imply 

J y = £daT4a)%±. (32) 

The above relation is a generalization of Eq. (J2J to non isothermal equilibria. 
The Ampere equation becomes 

dx 2 



da ^[-2£(A y (x) -a)+ 5] e -^ v(x )-ar+6 {Ay{ a»-a) _ (33) 

We can deduce a first order equation 

C+ r dak(a)e-^ a){A y- a)2+5( - A ^ x) - a) . (34) 



dx 

The vector potential is the solution of 

rM*) dA , . 

x = s (35) 

JMo) + jj dak(a) exp [-£(a)(A - a) 2 + 5 (a) (A - a)} ' 

The sign of the magnetic field is given by s = ±1, and C is an integration 
constant. The parameters k, £, and 5 are arbitrary functions of a. The 
solution has a physical meaning if it is defined for any value of x: the integral 
in Eq. (J33j) considered as a function of A y must be able to vary from — 00 to 
+00. 

Equation (pMjl can be modified in order to provide an integro-differential 
equation in B z (x). Let y defined by A y (y) = a be the new integration 
variable. Then, 

Bz (x) 2 = C+ T ^ J B 2 ( 2/ )My)e^ {s/)[ ^ Bz{u)du]2+5(2 ' ) ^ Bz(u) ' iul . (36) 

The same process applied on Eq. (J22|l provides 

»(*) = n + T dy^^^^e^C*^'^^ 3 '^. (37) 

These equations will be used in section III.B to analyse the asymptotic be- 
haviour of the solutions. 
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III. EXAMPLES AND PROPERTIES OF TAN- 
GENTIAL EQUILIBRIA 

A. Examples of elementary solutions with parameters 
relevant to space plasmas 

The trivial case of a uniform plasma can be easily recovered with several sets 
of parameters (including n g ^ and £ 7^ 0). The case of the Harris current 
sheet pHj, reviewed in details in ^2], can be recovered with £ = 0, v ^ and 
n(a) = n g So, where So is the Dirac distribution. 

Here are a few examples of equilibria where v = 0, n g = n c 5o, and 
where the temperature functions are constant. They correspond to the case 
g(p y ) = n c exp {—i]{p y /m) 2 ) already mentionned by Channell jTTj. The figure 
given in Channell's paper to illustrate this case displays the equilibrium of 
an evanescent plasma. We briefly show in this section a few examples of 
equilibria of finite size current layers in non evanscent plasmas, in order to 
emphasize a few interesting properties not discussed in Channell's paper. 

Figure 1 shows an example of a structure that is computed with pa- 
rameters typical of high altitude auroral plasmas : asymptotic magnetic field 
VC = 3000 nT, n = 10cm" 3 , T { = T e = 10 eV, and r] e = 1.53 10~ ie (m/s)- 2 = 
0.99?7 C , very close to the critical value imposed by Eq. (f2HJ) . A high nega- 
tive value n c = — 8cm~ 3 allows for a deep density depletion. The abcsissa is 
normalized to the ion Larmor radius which is pi = 152 m. We can see on 
Fig. 1 that the size of the structure is about 4p^ ~ 600 m. Such a narrow 
and deep structure (the density in the cavity is only 20% of the density of 
the surrounding plasma) can be used in a first approximation for the density 
depletions encountered in the high altitude auroral zone. 

The solution displayed in Fig. 2 corresponds to the limiting case for 
1], obtained for the weakly magnetized interplanetary plasma encountered 
in the solar wind at one astronomical unit : y/C = 5 nT, n = 10cm -3 , 
Ti = T e = 100 eV, and r] e = 1.53 10~ 17 (m/s)~ 2 = 0.99ry c . The ion Larmor 
radius is 280 km. We can notice that for the same density depletion, the 
variation of the magnetic field is much higher than in the case of a higly 
magnetized plasma. We also notice that the total size of the structure hardly 
exceed one ion Larmor radius. 

As the size is of the order of the ion Larmor radius, such solutions cannot 
be described through the fluid theories. Of course, we can also build very 
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large density structures, if we choose a small value of 77 in order to get a 
small value of £. Such large structures can more simply be described by fluid 
models. Figure 3 shows the solution with 77 = 10~ 20 (m/s) -2 ~ 10 _3 ^ c . The 
structure is much larger but has the same amplitude. Comparing Fig. 2 
and Fig. 3, we can notice that, contrary to the case of solitons, there is no 
correlation between the size and the amplitude of the structures. 

B. Asymptotic behaviour of the solutions 

We can caracterize two kinds of non trivial elementary solutions : those with 
77 7^ 0, and those where 77 = (and therefore v 7^ 0, otherwise the solution 
is trivial). Those of the first kind correspond to localized solutions, they 
are centered around an abscissa x such as A y (x) ~ a. Therefore, when the 
distribution is the infinite sum of elementary functions indexed by a as in Eq. 
([24)1. the parameter a indicates where the density n a (x) is non negligible. Let 
us consider an equilibria where, for x ~ 00, the magnetic field has a finite 
value. Then, the potential vector goes to infinity. Taking a,\ ~ — 00 and 
a 2 ~ +00 means that the plasma distribution will be influenced by some of 
those elementary functions even when x ~ 00. On the contrary, equilibria 
defined with finite values of a\ and a 2 (and B(oo) 7^ 0) means that for x ~ 00, 
the equilibria will converge toward the trivial solution of a uniform plasma 
with the density n that appears in Eq. 

The case of a null magnetic field at x ~ 00 (less useful for space plasmas 
applications) is different. Let us consider the case of _B(+oo) = 0. For 
x ~ +00, the vector potential has a constant value. Let us call Aoo this value. 
Then, as A cannot exceed A^, the main integral in Eq. ()35|) is divergent for 
A y (x) = Aoo (otherwise the equilibria would not be defined for any value of 
x). As the elementary solutions parametrized by a are non negligible only 
for abscissas x such that A y (x) is close to a, it is not necessary to consider 
values or a 2 that exceed notably A^ : they would have no influence on the 
value of the integral in Eq. ()35|). In brief, for equilibria with B(oo) = 0, the 
range of values where the parameter a is significant is bounded, and a 2 do 
not have to be infinite. 

In the Harris model, the magnetic field reverses and has two opposite 
finite values for x ~ — 00 and x ~ +00. We shall see that there exist non 
elementary solutions where the magnetic field has finite values on both sides, 
and that these values have not necessarily the same absolute value. 

We consider the case ai = —00 and a 2 = +00. Let u~ = \im x ^_ 00 u be 
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the limit of any physical value u when x tends to — oo, and u + = linx^+oo u. 
The integral in Eq. ()36|). where y\ = — oo and y 2 = +00 can be cut in three 
parts : an integral II from y\ to x — A, an integral 12 from a; — A to a; + A, and 
an integral 13 from x + A to y 2 - We will consider large values of x (positive 
or negative), and large values of A. 

For A large enough, y has large negative values, 

II ~ r A dyB;k-e W (-C(B;) 2 (x-y) 2 -5-B;(x-y)) (38) 
J yi 

~ -— [exp6~B~A-exp5~B~yi] (39) 

This function tends to zero for y\ — > —00 if 5~B~ < 0, and diverges if 
S~B~ > 0. In the case 5~ = 0, 

II ~ r _A ^ J B;A;-exp(-r(5;) 2 (x-i/) 2 ) (40) 

tends toward zero. Similarily, lim^^.oo 13 = if 5 + B + > 0. With the same 
kind of technique, we can show that, for large negative values of x, 

I* ~ fe"^ BJ exp^£, (41) 

where <7 B - is the sign of B~ . Provided that 5 + B + > 0, and 5~B~ < 0, the 
asymptotic expansion of Eq. (j3T?j) shows that 

{B-f = C + k-[^a (B7) exp (42) 



Computing an asymptotic value of 12 for x — > +00 bring a similar relation 
for (B+). The same method used with Eq. (J37j) provides the asymptotic 
value of the contribution of each species to the particle density: 

nr (s+) 2 

n+=n + N+a {Bt J — exp . (43) 

The equation for n~ is analogous. As d x n(oo) = 0, we deduce = J y = 0. 
There is no current density carried by the plasma for x ~ 00. This result 
can also be found with an asymptotic expansion of Eq. ()31|) . Moreover, Eq. 
(13 1|) shows that the ion and electron contribution to the current density have 
the same sign, their mean velocities are therefore opposite. As there is no 
current for x ~ 00, for each species v + = v ~ = 0. The plasma does not flow 
perpendicularly to the magnetic field at x ~ 00. 
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C. Examples of solutions with and without field reversal 



Figure 4 shows an example of equilibrium where the magnetic field amplitude 
and the density are not the same for x ~ — oo and x ~ +00. We have set 
constant values r\ = 1.53 lCT 17 (m/s)~ 2 close to r] c , u — 0, and a function 
n g ( a ) = 640atanh(1000a). The temperature is 10 eV. With C = (5nT) 2 , 
the asymptotic magnetic field amplitude is 5nT. The ion Larmor radius is 
Pi ~ 0.2 10 6 m. The variable a is homogeneous to B z times x, therefore 
1000a ~ x/ pi is of the order of magnitude of the absissa x divided by the 
ion Larmor radius. The factor 640 that multiplies the inverse hyberbolic 
tangent function was set heuristically in order to get a variation of density 
of the order of 60%. The relatively low value of \/~C = 5 nT allows for high 
relative variations of the magnetic field amplitude. The asymptotic values 
of the magnetic field and of the density fit the analytical formulas given in 
section III.B. 

Figure 5 shows an example where 77 = 0. We have chosen a function 
n g (a) = 4200e- {1000a)2 . The temperature is 10 eV, and VC = 5 nT. The value 
v = 10 -10 was chosen, because, in the exponential, the factor va ~ vB z x is 
again of the order of x/pi. We can see on Fig. 5 that the sign of the magnetic 
field changes for x ~ —00 and x ~ +00. The reversal occur for a value of 
A y that cancels the square root in Eq. (|33j) . Although the magnetic field is 
reversed on both sides, this is not a Harris equilibrium. 

D. Example of a non isothermal plasma cavity in the 
Earth auroral zone 

Are the cavities built in the Earth auroral zone in an equilibrium state of the 
plasma ? 

Hilgers et al. made a careful analysis of the Langmuir probe current 
measurements onboard the Swedish Viking satellite in low plasma densities. 
They showed taken at an altitude of 7000 km, of a deep auroral 

cavity. The density reach less than one particle per cubic centimeter, that 
is less than 10 % of the surrounding plasma density. The electron plasma 
temperature in the cavity (although not precisely measured) is of the order 
of 1 keV, compared to 1 eV outside the cavity. The boundaries of the cavity 
are sharp, Hilgers et al. measured 1.4 km, of the order of a few ion Larmor 
radii. The magnetic field amplitude is B zQ = 6800 nT. 

Apart from the temperature gradient, all these features have been qual- 
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itatively reproduced in the elementary equilibria shown in Fig. 1. In or- 
der to reproduce a cavity with sharp temperature gradients, an elemen- 
tary solution cannot be used. Let pi be the Larmor radius outside the 
cavity, it corresponds to the cold palsma (1 eV), its value is 21 m. In- 
side the cavity, the plasma temperature is 1 keV, and the ion Larmor ra- 
dius is Pi t hot = 660m= 32/?j. We choose an electron temperature function 
T e (a) = 1 + 1000exp(-(a/ J B z0 32p i ) 2 ) that varies from 1 to 1000 eV on a 
scale Ax ~ 32pj ~ Pi,hot- The density is controled through the constant 
scalar no = 10 cm -3 , and the n g function, n g (a) = —n g0 exp (— (a/ B z0 32pi) 2 ). 
The scalar n g o is set (heuristically) to 50000 in order to have one particle per 
cubic centimeter inside the cavity. The v function is null. The r] function is 
chosen is order to get a sharp cavity, that is with a value close to the r\ c func- 
tion: r\{a) = 0,98i] c (a) have values of the order of 10~ 15 (m/s)~ 2 . Actually, 
the equilibria is not very sensitive on the value of %]. Taking a lower value 
like rj(a) = 0.7rj c (a) (and a slightly lower value of n g0 ) bring similar results. 

The result, displayed on Fig. 6, is in quantitative agreement with the 
prescriptions given by the observations. The size of the gradient is 70 pi = 
1400m, as measured by Hilgers et al. 

We can notice on figure 6 that the magnetic field amplitude variation AB 
in the cavity is about 0.5 nT. A similar equilibria (not displayed) set with 
a lower value of the ambient magnetic field (B z0 = 300 nT) bring a larger 
variation: AB = 5 nT. 

We conclude from this short study that deep plasma cavities can be equi- 
librium structures of the Earth auroral plasma. Therefore, they won't be 
destroyed immediately after the extinction of their cause. It is however in- 
teresting to know if this equilibrium is stable. The theoretical treatement of 
this question goes beyhond the scope of this paper. I will only give a few 
hints in favour of the stability of the auroral cavities. 

The Auroral zone of the Earth is a radio source. The power of the emis- 
sions, called Auroral Kilometric Radiation (AKR), can reach 10 MW. The 
Viking spacecraft has gone through the sources of AKR: the waves are emit- 
ted inside the cavities, where they take their free energy from the hot rarefied 
plasma 7 s . The waves are strongly refracted at the edges of the AKR source, 
and it is probable that a part of the AKR is guided inside the cavities [II] . 
Auroral kilometric radiations can be observed on time scales of 10 minutes. 
The cavities that contain AKR sources are expected to last at least for the 
same duration. Therefore, two situations are possible : (1) the cavities are 
the consequence of a phenomena that lasts for tens of minutes, or (2) they 
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are generated by more transient phenomena and they are stable. 

An another hint about the stability of hot plasma cavities comes from 
numerical simulations carried by Genot et al. JO]. The simulations start 
with an isothermal (a few eV) auroral cavity. The cavity is stable. An 
Alfven wave is added. The interaction of the Alfven wave with the cavity 
triggers strong electron acceleration and turbulence. When the free energy is 
completely removed from the Alfven wave by the accelerated electrons, the 
plasma inside the cavity is hot, and the surrounding plasma is cold. We are 
in the situation shown in Fig 6. The simulations show that the heated cavity, 
which is still very deep and sharp, is stable. 

Most of the models about the generation of auroral cavities invoke strong 
kinetic Alfven Waves (SKAW) ^S] [113 • observed onboard the Freja and 
FAST satellites [TZj, ^H]- The associated magnetic field perturbation AB is 
of the order of 50 nT. But the observations made onboard FAST show that the 
SKAWs are mainly observed in the cusp and the polar cap boundary layer, 
while the cavities are observed in the lower latitude auroral region where the 
magnetic fluctuations AB do not exceed 5 nT. It is possible that the cavities 
are built in the polar cap boundary layer by SKAWs, they subsist after 
the disapearance of the Alfvenic perturbations and finally they are gently 
convected to lower latitudes by the large scale convection electric field. This 
suppose of course the stability of the auroral cavities. An other scenario is 
that the cavities are created through Field Line Resonance (FLR) JTUj directly 
in the auroral region. The AB ~ 5 nT (on a time scale that corresponds to 
the crossing of a cavity by a spacecraft) associated to the FLR is compatible 
with the value shown in Fig. 6. The equilibria shown in Fig. 6 may be a 
relevant short scale Vlasov- Maxwell description of the part of the FLR that 
is in the high altitude auroral zone (the FLR has been modeled on a global 
scale in the multifluid approximation up to now). 

IV. CONCLUSION AND FURTHER DEVLOPP- 
MENTS 

This paper presents a large class of monodimensional analytic kinetic equi- 
libria. They are based on particle distribution functions Eq. (jSJ) that depend 
on a set of almost arbitrary functions. These distribution functions are so- 
lutions of the Vlasov equation, and the vector potential can be computed 
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through the evaluation of the integral function given in Eq. (|35|) . In some 
particular cases, like the uniform plasma and the Harris current sheet, this 
integral is a combination of elementary mathematical functions. In the other 
cases, its numerical evaluation is straightforward. 

Unlike most of the kinetic tangential equilibria given in the litterature, the 
freedom in the choice of the density profile, altough more difficult to control, 
is almost as large as with bi-fluid models. The kinetic equilibria have the 
great advantage of giving a complete description of the distribution function 
(not provided with the fluid theory), and can describe equilibria where strong 
gradients develop on the scale of a few ion Larmor radii. These solutions 
can be used as initial conditions in particle in cell and Vlasov numerical 
simulations. 

Most of the authors who build models and are cited in this paper have 
based their elementary solutions on distribution functions that depend on the 
invariants of the motion. Some of them combined a small number of such 
elementary solutions and asserted that each of these elementary solutions 
correspond to a family of trapped particles. The distribution function given 
in Eq. (j3J) is the superposition, not of a finite number, but of a continuum 
of families of trapped particles. This is why we have a very large degree of 
freedom in the choice of the density and temperature profiles. For the first 
time, we show that such a supperposition of families of trapped particles lead 
to an analytically integrable equation, whose solution is given in Eq. (jHB} 

The equilibria discussed in section III.B where the plasma is uniform for 
x ~ oo can be compared to the solutions of the jump equations developed in 
the frame of the MHD theory. The equilibria developed in the present paper 
are characterized by jumps of the density and of the magnetic field, the 
normal velocity is null as well as the normal magnetic field. They therefore 
belong to the family of tangential discontinuities. But they concern only a 
sub-category: the magnetic field direction is uniform and the solutions have 
no velocity shear: the velocities are equal to zero at x ~ — oo and x ~ +00, 
even if they can take other values at finite distances. These restrictions are 
not required for general tangential discontinuities. 

Building a kinetic model, inspired from the present results, of a tangential 
discontinuity where the magnetic field can turn is straighforward and will be 
presented in a forthcomming paper in order to analyse experimental data 
provided by the Cluster satellites. 

A kinetic model of a tangential equilibria whith a velocity shift (x ~ — 00 
and x ~ +00) is not compatible with exact charge neutrality. A quasi neutral 
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model would do. It could be computed, as a first order perturbation in rii — n e 
added to the charge neutral equilibria presented in this paper. 

The examples given in the present paper come from the magnetospheric 
and solar wind physics because these media offers the opportunity of in situ 
observations of non collisional astrophysical plasmas. Tangential equilibria 
can be used to describe the magnetopause (present around all the magnetized 
planets with an atmosphere), the distant neutral sheet, the Earth auroral 
cavities, and some density fluctuations of the solar wind [2U| . 

But tangential equilibria do not exist only in the terrestrial environnment. 
They exist whenever plasmas with different origins and velocities meet. Such 
situations exist in many astrophysical plasmas. Theoretical models (based 
on MHD and plasma multifluid theories) of the boundary of the Heliosphere 
show that it is constituted of a termination shock, followed by the heliopause, 
that is a tangential discontinuity, an hydrogen wall, and possibly an helio- 
spheric shock [21]. Tangential discontinuities also exist around other stars. 
Miiller et al. [221 studied the interaction of the very active binary star A An- 
dromedae, and of the neaby star e Indi with the interstellar medium. They 
found, as in the case of the heliopause, the existence of four boundaries, one 
of them, the asteropause being a tangential discontinuity. It is clear that 
tangential layers exist in other astrophysical plasmas where they may play 
a very important role, as frontiers, or in acceleration, heating and radiative 
processes. 

The studies of the magnetospheric plasmas have shown that most of the 
heating and acceleration phenomena can be only explained through kinetic 
processes. So far, most of the remote astrophysical plasmas have been stud- 
ied through MHD or multifluid theories. The understanding of acceleration 
processes in remote astrophysical plasma might as well require kinetic mod- 
els. The equilibria presented in this paper might be a good start for some of 
those studies. 
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Figure 1: FIG. 1. A deep and narrow density depletion in a highly magne- 
tized plasma. The parameters of this equilibrium are given in Sec. III. A. 
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Figure 2: FIG. 2. A deep and narrow density depletion plasma with a weak 
magnetic field. The parameters are the same as for Fig. 1, exepct for C that 
is much weaker. 21 
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Figure 3: FIG. 3. A deep and large density depletion, the only difference 
with Fig. 2 is a larger value of rj. See Sec. III. A.. 
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Figure 4: FIG. 4. An example of non symetric equilibrium where n g is a 
function of a. See Sec. III.C for details. 
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Figure 5: FIG. 5. An exemple of an equilibrium with r] = 
of the magnetic field. Details are given in Sec III.C. 
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and a reversal 
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Figure 6: FIG. 6. An example of a deep plasma cavity containing hot elec- 
trons (1 keV) surrounded by a cold (1 eV) highly magnetized plasma. See 
section III.D for details. 25 



